Back

Crop Science

Wiley

Preprints posted in the last 90 days, ranked by how well they match Crop Science's content profile, based on 18 papers previously published here. The average preprint has a 0.02% match score for this journal, so anything above that is already an above-average fit.

1
Single-kernel near-infrared spectroscopy enables haploid kernel sorting in field and sweet corn using high-oil haploid inducers across diverse donor-inducer combinations

Sharma, S.; Gustin, J. L.; Frei, U. K.; Settles, A. M.; Lübberstedt, T.; Resende, M. F. R.; Hershberger, J.

2026-07-24 plant biology 10.64898/2026.07.23.740370 medRxiv
Top 0.1%
13.1%
Show abstract

Key messageA single-kernel near-infrared reflectance spectroscopy-based sorter can effectively identify haploid kernels for doubled haploid production in field and sweet corn backgrounds. Doubled haploid (DH) technology significantly shortens the breeding cycle for developing homozygous inbred lines in maize (Zea mays). Manual sorting of haploids from a larger bulk of hybrid kernels in an induction cross is a major bottleneck in DH development. Automated systems based on near-infrared (NIR) reflectance spectroscopy can be valuable tools for rapid haploid sorting, provided that sorting accuracy is sufficient for incorporation into the DH process. In this study, we evaluated the accuracy of a custom-built single-kernel NIR (skNIR) sorter for classifying haploid kernels from 12 high-oil haploid induction populations generated from two sweet corn and two field corn donors and four high-oil haploid inducers (HOHIs). We evaluated several general classification models that can be applied without population-specific recalibration or prior genotyping, including models that classified haploids based solely on predicted oil content, as well as multivariate methods that used all wavelengths of the NIR spectra. The highest classification accuracy was obtained using a general multivariate support vector machine (SVM) model. When combined with the two best-performing HOHIs, the general SVM model accurately sorted induction populations from two of the three donor backgrounds crossed with these inducers. Two oil-based methods showed less accurate classification than the multivariate SVM model, due to overlapping oil content distributions across the two kernel classes. Overall, this study demonstrates effective skNIR-based sorting of haploid kernels from diverse induction populations using a single general model. The practical deployment of this instrument in maize breeding programs is discussed.

2
Long-term realized genetic gain and population dynamics under genomic selection in Brazilian cassava germplasm

de Freitas, G. M.; Certuche, D. C. S.; Jannink, J.-L.; De Oliveira, E. J.; Garcia, A. A. F.

2026-07-22 genetics 10.64898/2026.07.18.739356 medRxiv
Top 0.1%
11.2%
Show abstract

Genomic selection has become an important strategy in cassava breeding, enabling faster selection cycles and sustained genetic progress. Despite its widespread adoption, long-term evaluations integrating predictive performance, realized genetic gain, and genetic diversity remain scarce, particularly in clonally propagated crops. We present a comprehensive assessment of genomic selection outcomes in the Brazilian cassava breeding program across four recurrent selection cycles (C0 to C3) implemented between 2011 and 2024, using historical phenotypic and genomic data from 210 multi-environment trials. Predictive ability of genomic best linear unbiased prediction models ranged from low to moderate, depending on the traits genetic architecture and heritability. Prediction accuracies were highest in early cycles (C0 and C1) and showed modest declines in later cycles (C2 and C3). Root yield, shoot yield, plant height, starch content, and dry matter content exhibited stable predictive performance across cycles, with a gradual reduction in RMSE, indicating improved model calibration as training populations expanded. Regression analyses of genomic estimated breeding values revealed significant realized genetic gains for most yield-related traits. In contrast, dry matter content and starch content exhibited small, non-significant negative trends, consistent with known unfavorable genetic correlations with yield. Targeted reductions in plant architecture scores reflected deliberate selection for ideotypes suited to mechanized production systems. At the same time, analyses of genetic diversity revealed a slight decrease in observed heterozygosity, with higher values in the most advanced selection cycle. These results provide an integrated framework for monitoring predictive performance, realized genetic gain, and population genetic dynamics under long-term genomic selection. Collectively, they offer valuable insights into balancing short-term genetic improvement with long-term sustainability and support the development of strategies to optimize selection decisions, breeding planning, and population management in Brazilian cassava breeding programs.

3
From long-term progeny trials to genomic selection: empirical prediction and simulation-guided redesign of Scots pine breeding in Germany

Degen, B.; Leite Montalvao, A. P.; Schneck, V.

2026-07-26 genetics 10.64898/2026.07.21.739892 medRxiv
Top 0.1%
6.7%
Show abstract

Long breeding cycles constrain genetic gain in Scots pine, while mature progeny trials can provide reference populations for genomic selection. We combined phenotypic records from two German trials established in 1990 with dense SNP data. We complemented these data with an offspring-level parentage audit, duplicate filtering, and simulations of breeding strategies. Of 1,986 phenotype-matched genotyped trees, 1,668 formed the parentage-green set and 1,651 remained after exclusion of 17 near-duplicate samples. The final population comprised 1,521 supported controlled-cross offspring and 130 mother-known, father-unknown offspring from 87 progeny labels. Across 14 trait-by-age measurements, GBLUP gave the clearest improvements for diameter and volume from age 20 onwards, whereas height was more mixed. At age 35, ABLUP versus GBLUP heritability was 0.155 versus 0.235 for diameter, 0.265 versus 0.263 for height, and 0.172 versus 0.236 for volume. Leave-progeny-out predictive ability increased from 0.205 to 0.238, 0.249 to 0.265, and 0.203 to 0.239, respectively. SNPscan_breeder simulations compared phenotypic selection, progeny testing, cross-generation genomic selection, and genomic selection with phenotypic thinning under five diversity variants. Progeny testing produced the greatest cumulative gain, but its 45-year cycle reduced annual response. Under the base assumptions, genomic selection with phenotypic thinning gave the highest annual gains for height and fungal resistance, whereas pure genomic selection gave the highest annual diameter gain. Genomic strategies accumulated more kinship than conventional strategies, although a {lambda} = 0.10 kinship penalty improved founder retention with little loss of gain; no mitigation option simultaneously maximised gain, prediction accuracy, and diversity. These results support genomic shortlisting within a field-tested programme with explicit reference-population updating and diversity management.

4
Identification of heterotic group-specific haplotypes and impact of residual inbreeding on grain yield of maize elite hybrids

Kadoumi, R.; Heslot, N.; Henriot, F.; Murigneux, A.; Berton, M.; Moreau, L.; Charcosset, A.

2026-06-21 genetics 10.64898/2026.06.15.732226 medRxiv
Top 0.1%
6.2%
Show abstract

Modern hybrid maize (Zea mays L.) breeding programs are based on the management of distinct complementary heterotic groups to maximize heterosis in high-performing hybrids. This practice lowers shared genetic segments and increases divergence between groups to limit inbreeding in hybrids. However, most breeding programs have not always enforced strict separation between heterotic groups in the past. Competitor commercial hybrids were notably a common elite germplasm source for inbred development, which would diminish divergence between groups. This study proposes a new haplotype-based approach to assess hybrids residual inbreeding based on parental similarity. The new haplotype method has a stronger significant negative effect on hybrids grain yield than raw SNP data. Evaluation of modern experimental hybrids uncovered related inbreds contributing to superior rates of residual inbreeding. Analysis of these inbreds revealed haplotype transfers between heterotic groups, originating notably from the use of a Stiff Stalk-Iodent commercial hybrid as breeding starts material in both Stiff Stalk and Non-Stiff Stalk breeding populations. The introduction of this intergroup parent generated heterotic-group-specific haplotype migration between crossing pools. These fragments caused significant genome-wide residual inbreeding in experimental hybrids across selection cycles. This study highlights the necessity for accurate evaluation of external sources of diversity to minimize haplotype transfers and admixture between crossing pools. We demonstrate the consequences of using commercial hybrids in inbred development, particularly regarding residual inbreeding, and their effects on hybrid performance. Insights from these results can assist breeders in optimizing the choice of parents for introducing genetic diversity in a reciprocal recurrent selection scheme. KEY MESSAGEHaplotype-based hybrids parental similarity better predicts grain yield than marker-based identity-by-state. Utilization of commercial hybrids as breeding start material resulted in higher hybrid residual inbreeding even after several selection cycles

5
Comparison of localGEBV and Optimal Haplotype Stacking Fitness Functions using a Novel R Package: HapSelect

Shaffer, W.; Papin, V.; Carter, Z.; Brunner, S. M.; Tong, J.; Villiers, K.; Robinson, H.; Voss-Fels, K.; Hayes, B. J.; Hickey, L.; Dinglasan, E.

2026-07-13 genetics 10.64898/2026.07.08.737160 medRxiv
Top 0.1%
5.4%
Show abstract

Haplotype-based breeding strategies have emerged as promising approaches to maximize long-term genetic gain by identifying complementary parental combinations while maintaining genetic diversity. However, these methods typically require phased genotypes and more intensive workflow pipelines and skillsets. We developed a novel local genomic estimated breeding value (localGEBV) fitness function with similar intent to the optimal haplotype stacking (OHS) framework fitness function and implemented both in the novel R package, HapSelect. Our aim was to evaluate whether phased haplotypes provide additional benefit over the more easily available dosage-based unphased genotypes in highly inbred crops. A subset of bread wheat nested association mapping (NAM) population comprising 444 lines genotyped with 6,054 DArT-Seq markers was analysed. Marker effects were estimated using rrBLUP, localGEBV and haplotype effects were calculated across linkage disequilibrium-defined haploblocks, and genetic algorithms (GA) were used to identify optimal sets of 30 founders using either a localGEBV derived fitness function with unphased, dosage inputs or the OHS fitness function with phased inputs. Selected parental sets were compared with conventional truncation selection (TS) through 150 generations of forward simulation. The OHS fitness function achieved a marginally greater optimized ultimate GEBV than the localGEBV fitness function during GA optimization, with only 18 of the 30 selected founders overlapped between the two methods. Despite these differences, forward simulations demonstrated nearly identical long-term genetic gain for localGEBV and OHS-selected founders, with both approaches outperforming conventional truncation selection by maintaining greater genetic diversity and delaying the genetic plateau. The minimal difference between localGEBV and OHS is likely attributable to the high homozygosity of the population, where localGEBV and haplotype effects are nearly confounded. These results demonstrate that dosage-based localGEBV provides a practical alternative to phased haplotype approaches for parent selection in inbred crops, substantially simplifying genomic workflows while maintaining long-term breeding performance. Future work should evaluate these methods in more diverse inbred populations and outbred species, where great haplotypic diversity may increase the advantage of true haplotype-based optimizations.

6
Cross Potential Selection for Multiple Traits Considering the Progeny Distribution of Future Inbred Lines in Plant Breeding Programs

Sakurai, K.; Moreau, L.; Mary-Huard, T.; Charcosset, A.; Iwata, H.

2026-06-08 genetics 10.64898/2026.06.02.729654 medRxiv
Top 0.1%
5.3%
Show abstract

In plant breeding, it is often necessary to improve a target trait while maintaining other essential traits within desirable ranges. When genetic relationships exist among these traits, improvements in the target trait may lead to undesirable changes in essential traits, complicating cross selections. In such cases, it is critical to select cross-pairs that are expected to produce progeny that satisfy the requirements for all traits. The progeny distribution of each crossing pair can be predicted using the estimated genotypic values and genetic (co)variances of the target and essential traits. By utilizing this distribution, the probability of generating progeny that satisfy predefined trait requirements can be evaluated, allowing a direct comparison of alternative crosses. In this study, we developed Cross Potential Selection for Multiple Traits (CPS-MT), a breeding strategy designed to improve a target trait while maintaining one or more essential traits within desirable ranges. CPS-MT extends the original Cross Potential Selection (CPS) framework to explicitly handle trade-offs between traits under genetic correlations. We evaluated the performance of CPS-MT through simulations involving four types of genetic relationships and two genetic causal factors between traits, resulting in seven scenarios. Across all scenarios, CPS-MT consistently improved the likelihood of obtaining desirable progeny, indicating that CPS-MT provides a practical and effective framework for cross selection under multi-trait constraints in breeding programs. Article SummaryThis study developed Cross Potential Selection for Multiple Traits (CPS-MT), a new breeding strategy designed to improve a target trait while maintaining one or more essential traits within desirable ranges. CPS-MT evaluates crossing pairs by predicting progeny distributions based on estimated genotypic values and genetic covariances, enabling direct comparison of alternative crosses under multi-trait constraints. Through simulations incorporating four types of genetic relationships and two causal factors (seven scenarios), CPS-MT consistently increased the likelihood of obtaining progeny that satisfied the predefined trait requirement. These results indicate that CPS-MT provides a practical, robust framework for target trait improvement under trait constraints.

7
Assessing the Role of Marker Density and Minor Allele Frequency on Machine Learning Driven Genomic Selection Accuracy in Grapevine

Francisco, F. R.; de Oliveira, G. L.; Niederauer, G. F.; Fritsche-Neto, R.; Souza, A. P. d.; Furlan, M. F. M.

2026-07-17 genetics 10.64898/2026.07.11.737951 medRxiv
Top 0.1%
5.0%
Show abstract

Although grapevine (Vitis spp.) is among the oldest and most economically significant fruit species globally, its genetic improvement faces major bottlenecks due to long juvenile periods and extended cycles for phenotypic evaluation. In this context, genomic selection (GS) has emerged as an effective alternative to traditional selection, offering a robust framework to optimize breeding programs by significantly reducing generation intervals while enhancing predictive accuracy (PA) in early generations and expected genetic gains (EGGs). Nevertheless, factors such as minor allele frequency (MAF) and population size can significantly affect predictive models, even to the point of making their use unfeasible in breeding programs. In this context, this study evaluated the effect of data dimensionality reduction on GS accuracy by selecting single-nucleotide polymorphisms (SNPs) based on MAF thresholds. The experimental design tested the predictive capacities of four machine learning (ML) algorithms (ElasticNet, K-Neighbors, Support Vector Machine Regression, and XGBoost) alongside the conventional Genomic Best Linear Unbiased Prediction (gBLUP) model. These were validated using three SNP datasets (11,115, 9,494, and 6,100 markers) filtered by MAF levels of 0.05, 0.1, and 0.2 across six genetic traits, and EGGs were compared between conventional breeding and GS via the breeders equation. The results revealed that the ML models exhibited remarkable stability, with no significant differences in PA across the different MAF-based SNP densities, except for berry length, which showed a substantial difference with XGBoost at an MAF of 0.2. Conversely, gBLUP demonstrated high sensitivity to dimensionality reduction, with its performance significantly impacted by MAF filtering across all the traits. These results suggest that compared with traditional GS models that rely on a genomic kinship matrix, ML-based approaches offer greater flexibility in feature reduction. Additionally, compared with chemical traits, morphological traits generally had greater predictive ability. Furthermore, every GS model provided estimated genetic gains superior to traditional breeding, with improvements ranging from an 8.90-fold increase in berry length to a 2.86-fold increase in total soluble solids, confirming that GS integration is promising for enhancing breeding efficiency in grapevines.

8
Impact of Reduced Chlorophyll Levels in Leaves on Soybean Yield, Seed Composition, Pod/Seed Photosynthesis, and Chlorophyll Levels in Pod and Seed Tissues

Jones, S. I.; Stutz, S. S.; Atalay, E.; Wang, Y.; Ort, D. R.; Cho, Y. B.

2026-08-19 plant biology 10.64898/2026.08.14.744892 medRxiv
Top 0.1%
4.7%
Show abstract

Soybean, a widely cultivated leguminous crop valued for its protein, amino acids, and oil, faces the challenge of maintaining protein levels, which have an inverse correlation with yield. Reducing leaf chlorophyll levels could increase seed protein levels without compromising yield; however, this is yet to be tested. Therefore, to understand the impacts of low chlorophyll mutations on soybean yield and seed composition, we screened and compared 25 low chlorophyll soybean mutants to their 11 dark green parents. PI548210 (Lincoln mutant) demonstrates a higher concentration of protein without affecting yield compared to its dark green parent PI548362 (Lincoln), suggesting it as a good candidate for further large-scale field trials. PI547555 (Y11/y11, Clark mutant) demonstrates a lower concentration of oil without impacting yield, alongside lower gross photosynthesis, but with chlorophyll levels in the pod and seed tissues that are comparable to its dark green parent PI548533 (Clark). These findings are consistent with the oil concentration of the soybean being influenced by pod and seed photosynthesis, which is correlated with pod height and row spacing. Chlorophyll levels in the leaf do not necessarily correlate with those in the pod and seed of low chlorophyll mutants, possibly due to substantially lower expression of chlorophyll synthesis genes in the pod and seed. SIGNIFICANCEO_LIPI548210 (Lincoln mutant), one of twenty-five low chlorophyll soybean mutants, demonstrates a higher concentration of soybean protein without affecting yield compared to its dark green parent (Figure 1 and Table 1). C_LIO_LIPI547555 (Y11/y11, Clark mutant), a low chlorophyll soybean mutant, demonstrates a reduced concentration of soybean oil without impacting yield, alongside lower gross photosynthesis in pod and seed tissues compared to its dark green parent (Figures 3 and Table 2). These findings suggest that the oil concentration of the soybean is influenced by pod and seed photosynthesis, which is in turn influenced by pod height and row spacing (Figure 2). C_LIO_LIChlorophyll levels in the leaf do not necessarily correlate with those in the pod and seed of low chlorophyll mutants, possibly due to substantially lower expression of chlorophyll synthesis genes in the pod and seed (Figure 5-6). C_LI O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=84 SRC="FIGDIR/small/744892v1_fig1.gif" ALT="Figure 1"> View larger version (55K): org.highwire.dtl.DTLVardef@4282dcorg.highwire.dtl.DTLVardef@9d565forg.highwire.dtl.DTLVardef@1918292org.highwire.dtl.DTLVardef@1359b1_HPS_FORMAT_FIGEXP M_FIG O_FLOATNOFigure 1.C_FLOATNO Two low chlorophyll mutants are as healthy as their dark green parents. Lincoln and its low chlorophyll mutant, left; Clark and its low chlorophyll mutant, known as Y11/y11, right. It can be seen by eye that the plants have low chlorophyll (light green/yellow leaves) but a similar growth habit to their dark green parents. See Supplemental Figures 1-4 for contrast, where low chlorophyll mutants are stunted in growth compared to their dark green parents. C_FIG O_TBL View this table: org.highwire.dtl.DTLVardef@657ec9org.highwire.dtl.DTLVardef@166e75borg.highwire.dtl.DTLVardef@df23c7org.highwire.dtl.DTLVardef@1a60124org.highwire.dtl.DTLVardef@194ed96_HPS_FORMAT_FIGEXP M_TBL O_FLOATNOTable 1.C_FLOATNO O_TABLECAPTIONComparison of seed yield, weight, seed composition between low chlorophyll mutants and their dark green parents. ANOVA is used with linear mixed model (random effect = block, fixed effect = variety). Least squares mean is used to compare. For yield and seed composition, N=4 blocks. For leaf chlorophyll (SPAD), N=40. Yield is average yield per plant (g). n.s. = not significant. C_TABLECAPTION C_TBL O_FIG O_LINKSMALLFIG WIDTH=179 HEIGHT=200 SRC="FIGDIR/small/744892v1_fig3.gif" ALT="Figure 3"> View larger version (26K): org.highwire.dtl.DTLVardef@7a368aorg.highwire.dtl.DTLVardef@192b8f0org.highwire.dtl.DTLVardef@1abb738org.highwire.dtl.DTLVardef@89e978_HPS_FORMAT_FIGEXP M_FIG O_FLOATNOFigure 3.C_FLOATNO Light response curve of low chlorophyll mutant (Y11/y11, PI547555) and its parent (Clark, PI548533). Rates of net and gross photosynthesis of low chlorophyll (white) and dark green parents (black) pods under field conditions. Each dot represents a value (n=4) {+/-}SE. We assumed that the seeds greatly inhibited the transmittance of light through the pod and used photosynthetic photon flux density for a single-side. C_FIG O_TBL View this table: org.highwire.dtl.DTLVardef@3f0528org.highwire.dtl.DTLVardef@16ba712org.highwire.dtl.DTLVardef@a5ab2aorg.highwire.dtl.DTLVardef@889254org.highwire.dtl.DTLVardef@3efa4f_HPS_FORMAT_FIGEXP M_TBL O_FLOATNOTable 2.C_FLOATNO O_TABLECAPTIONPod photosynthetic parameters for low chlorophyll mutant (Y11/y11, PI547555) and its parent (Clark, PI548533). Photosynthesis was measured 1 September through 15 September 2021 at the University of Illinois Energy Farm in Urbana, IL, USA. The statistical analysis was done using ANOVA with linear mixed model (alpha=0.05). N=4 {+/-} SEM for Clark and N=3 {+/-} SEM for Y11. C_TABLECAPTION C_TBL O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=130 SRC="FIGDIR/small/744892v1_fig2.gif" ALT="Figure 2"> View larger version (23K): org.highwire.dtl.DTLVardef@a36c26org.highwire.dtl.DTLVardef@1116c8forg.highwire.dtl.DTLVardef@ee5e61org.highwire.dtl.DTLVardef@1766712_HPS_FORMAT_FIGEXP M_FIG O_FLOATNOFigure 2.C_FLOATNO Low chlorophyll mutant (Y11/y11, PI547555) and its parent (Clark, PI548533) differ in concentration of seed oil, which interacts with height of pod and row spacing. The box plots show the median (central line), the lower and upper quartiles (box) and the minimum and maximum values (whiskers). The statistical analysis was done using ANOVA with linear mixed model (n=3 blocks, alpha=0.05). Least squares mean is used to compare. N.s., non- significant in the analysis. A. Concentration of oil in low chlorophyll mutant seeds from the upper canopy decreased by 4% compared to the dark green parent (18.2% vs 19%) while there was no difference between them in the seeds from the lower canopy (20.2% vs 20.6%). B. Schematic layout of 2013 field setting showing two different row spacings. C. Concentration of oil in low chlorophyll mutant decreased by 2% in 38cm spacing (21.4% vs 22%) while there was no difference in 19cm spacing (21.3% vs 21.7%) in 2013 field. C_FIG O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=162 SRC="FIGDIR/small/744892v1_fig5.gif" ALT="Figure 5"> View larger version (22K): org.highwire.dtl.DTLVardef@68e508org.highwire.dtl.DTLVardef@94a6ccorg.highwire.dtl.DTLVardef@152a187org.highwire.dtl.DTLVardef@1eae137_HPS_FORMAT_FIGEXP M_FIG O_FLOATNOFigure 5C_FLOATNO (greenhouse). Correlation between the level of leaf chlorophyll (x-axis: SPAD reading) and the level of immature pod or seed chlorophyll (y-axis, mg/g DW). Line represents the linear regression model. R-squared is a coefficient of determination, the percentage of the response variable variation that is explained by the linear model. Pod is labeled by the fresh weight of seeds it contained. A. Level of chlorophyll of 25-100mg pod (n=18). B. Level of chlorophyll of 100-200mg pod (n=17) . C. Level of chlorophyll of 25-100mg seed (n=17). D. Level of chlorophyll of 100-200mg seed (n=20). C_FIG O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=180 SRC="FIGDIR/small/744892v1_fig6.gif" ALT="Figure 6"> View larger version (28K): org.highwire.dtl.DTLVardef@167fd88org.highwire.dtl.DTLVardef@361472org.highwire.dtl.DTLVardef@786325org.highwire.dtl.DTLVardef@1b53855_HPS_FORMAT_FIGEXP M_FIG O_FLOATNOFigure 6.C_FLOATNO Levels of gene expression in chlorophyll synthesis pathway. A. CHL common pathway genes; Glutamyl-tRNA reductase (GluTR). Glutamate 1- semialdehyde aminotransferase (GSA-AT). ALA dehydratase (ALAD). Uroporphyrinogen III synthase (UROS). Uroporphyrinogen III decarboxylase (UROD). Protoporphyrinogen IX oxidase (PPO). B. Mg branch; Mg-chelatase (Mgch). Magnesium-protoporphyrin IX monomethyl ester cyclase (MPEC). Protochlorophyllide reductase (POR). 3,8-divinyl protochlorophyllide a 8-vinyl-reductase (4VCR). Heme pathway; Ferrochelatase (FECH). Heme oxygenase (HO). Phytochromobilin synthase (HY). Data come from Severin et al (2010). RPKM, reads per kilobase per million mapped reads. DAF, days after flowering. The source seed is experimental line A81-356022 which was generated by introgressing G. soja (PI468916) into G. max (A81-356022). C_FIG

9
Enhancing predictive accuracy of yield traits in cassava through multi-trait genomic prediction

de Freitas, G. M.; Certuche, D. S.; Jannink, J.-L.; de Oliveira, E. J.; Garcia, A. A. F.

2026-07-06 genetics 10.64898/2026.07.01.735838 medRxiv
Top 0.1%
4.4%
Show abstract

Multi-trait genomic prediction offers a practical route to improve selection for costly, complex traits in clonally propagated crops such as cassava. In a Brazilian breeding panel of 1,078 cassava clones genotyped with 25,923 SNPs and phenotyped for six agronomic traits, we compared single-trait (ST) and multi-trait (MT) GBLUP models. Stage-wise mixed models produced BLUEs that fed into ST and MT-GBLUP. We tested five cross-validation schemes that mimic breeder realities: ST baseline (CV1); naive all-traits MT prediction for unphenotyped candidates (CV2); MT prediction using auxiliary trait phenotypes in the test set (CV3); and two sparse-phenotyping regimes with missingness by trait (CV4) or by clone (CV5) at 25%, 50%, and 75% levels. The main results were that, under the ST baseline (CV1), predictive ability ranged from 0.50 for DMC and 0.45 for FRY down to 0.13 for Le.Dis. A naive full MT model (CV2) performed approximately on par with ST-GBLUP. In contrast, MT designs (CV3) that included informative auxiliary traits, such as shoot yield and combinations with plant vigor and leaf disease severity, yielded small gains for DMC with predictive ability of approximately 0.51 (+2%), while FRY predictive ability increased to approximately 0.65 (+44%), accompanied by RMSE reductions for FRY up to approximately 13.5% (e.g. RMSE approximately 6.2). Sparse-phenotyping simulations (CV4/CV5) demonstrated that MT models sustain or even improve predictive ability under realistic missing-data regimes (PA {approx} 0.62 - 0.65). Selection concordance between MT and ST top-10% sets was generally high (>0.80), and MT configurations produced measurable improvements in expected selection response and genetic gain per cycle for several target traits. These results indicate that strategically implemented MT-GBLUP, using a small set of biologically and operationally informative auxiliary traits and optimized sparse phenotyping, can materially increase predictive accuracy and selection efciency for economically critical cassava traits while reducing phenotyping burden.

10
Novel quantitative trait loci conferring broad-based resistance to root-knot nematodes in lima bean (Phaseolus lunatus)

Tajima, A. M.; Matthews, W. C.; Duong, T.; Khanh, T. D.; Baniya, A.; Penmetsa, R. V.; Parker, T.; Farmer, A.; English, S.; Diepenbrock, C.; Gepts, P.; Roberts, P. A.; Huynh, B.-L.

2026-07-09 plant biology 10.64898/2026.06.30.735594 medRxiv
Top 0.1%
4.0%
Show abstract

Lima bean (Phaseolus lunatus) is a broadly adapted, economically important leguminous crop and a susceptible host of root-knot nematodes (Meloidogyne spp.; RKN), which are a devastating plant pathogen in agricultural systems worldwide. To date, there have been few studies to elucidate the genetic determinants of RKN resistance in lima beans. Understanding the genetic mechanisms underlying resistance is essential for improving resistance traits and incorporating them into lima bean breeding programs. To assist in marker-assisted selection, we aimed to identify and map quantitative trait loci (QTLs) conferring RKN resistance-related traits. Three recombinant inbred line (RIL) populations were used in this study. Three populations were derived by crossing two RKN-resistant parents with the same RKN-susceptible parent and with each other. All populations were genotyped using genome-wide single-nucleotide polymorphism (SNP) markers. Each population was screened for root galling (RG) and RKN egg reproduction (ER) in response to M. incognita and M. javanica in greenhouse experiments. Three major QTLs were detected and mapped on chromosome Pl04 (QRk-pl04.1), Pl05 (QRk-pl05.1) and Pl10 (QRk-pl10.1) across populations. Among them, QRk-pl05.1 and QRk-pl10.1 affected levels of RG and ER of both RKN species, while QRk-pl04.1 suppressed root galling and reproduction responses of M. incognita but not of M. javanica. These chromosomal regions defined by flanking markers will help guide marker-assisted breeding and gene discovery for broad-based RKN resistance in lima beans.

11
Introducing PHJ Media: A Unique Machine Learning -Driven Basal Formulation to Overcome Recalcitrance for Multi-Genotype Micropropagation of Cannabis sativa L.

Pepe, M.; Hesami, M.; Jones, M.

2026-07-15 plant biology 10.64898/2026.07.14.738465 medRxiv
Top 0.1%
3.3%
Show abstract

Applications of tissue culture are critical for Cannabis sativa L. (cannabis), supporting clonal propagation, germplasm preservation, pathogen elimination, among other biotechnological applications. However, extensive genetic diversity associated with cannabis results in highly variable responses to in vitro conditioning, and no consensus basal media formulation exists to support reproducible micropropagation across genotypes. To address these limitations, a hybridized ensemble-NSGA-II approach was employed for concurrent optimization of individual media components to create a species specific, cultivar inclusive basal salt formulation for cannabis micropropagation. The resulting PHJ media represents a unique formulation that overcomes recalcitrance across a wide array of cannabis cultivars, facilitating improved growth and uniformity for the nine cultivars used in its development and validation. These results remain consistent from explant initiation through multiple rounds of subculture. The ability of PHJ to overcome genotypic recalcitrance is telling of its potential applicability with an array of plant species beyond cannabis. Additionally, robust performance both with and without plant growth regulators underscores the plausible use of PHJ for diverse applications beyond standard micropropagation. Ultimately, this cultivar-inclusive basal medium demonstrates utility for both scientific research and industrial-scale operations.

12
Genetic Diversity and Population Structure of Maize Doubled Haploid Lines from Drought and Low Nitrogen Tolerant Populations

Ehemba, G. L.; Ifie, B. E.; DAS, B.; Abu, P.; Adjei, E. A.; Ayenan, M. A. T.; Garcia-Oliveira, A.; Ribeiro, P.; Manilal, W.; Tongoona, P.; Danquah, E. Y.

2026-08-13 genetics 10.64898/2026.08.06.743182 medRxiv
Top 0.1%
3.3%
Show abstract

Understanding the genetic diversity and population structure of breeding materials is essential for developing stress-resilient cultivars. In tropical maize, where drought and low soil nitrogen (low N) severely limit productivity, continuous development of tolerant varieties remains a priority. This study assessed the genetic diversity and population structure of 250 doubled haploid lines (DHLs) derived from five drought- and low N-tolerant tropical populations. Genotyping was performed using mid-density DArTseq markers, yielding 3,305 high-quality SNPs for analysis. Results revealed a moderate level of diversity among the DHLs, with an average genetic distance of 0.39, a polymorphism information content (PIC) of 0.33, and a minor allele frequency (MAF) of 0.29. These values reflect substantial allelic variation, important for identifying complementary parental combinations in hybrid development. Discriminant analysis of principal components (DAPC) grouped the DHLs into five distinct clusters, largely corresponding to their source populations, although some admixture was observed. This indicates that while the genetic backgrounds of the source populations were mostly retained, recombination introduced useful variation. Overall, the clear population structure and high diversity observed among these DHLs provide a strong genetic foundation for future maize improvement. These lines represent valuable resources for heterotic group formation, hybrid development, and recurrent selection schemes aimed at enhancing drought and low nitrogen tolerance in tropical maize.

13
Far-red timing uncovers cultivar-dependent yield and bolting responses in vertical-farm spinach (Spinacia oleracea L.)

McGovern, C.; Adrio, M.; Aliki, H.; Vichos, R.; Powell, W.; Sharma, R.

2026-07-13 plant biology 10.64898/2026.07.10.737849 medRxiv
Top 0.1%
2.7%
Show abstract

Far-red light (FR; 700-750 nm) is increasingly incorporated into controlled-environment lighting because it can improve photosynthetic efficiency when combined with comparatively shorter wavelengths. In long-day leafy crops such as spinach, however, FR may also promote the transition from vegetative to reproductive growth and thereby reduce marketable yield. Most studies have evaluated FR fraction, intensity or end-of-day exposure, whereas the developmental timing of FR has rarely been tested, particularly in spinach. Here, we evaluated six commercial spinach cultivars (Amador, Harp, Renegade, Responder, Rubino and Santa Cruz) in an indoor vertical farm under a common red-green-blue background (PPFD 260-264 {micro}mol m-{superscript 2} s-{superscript 1}, 12 h photoperiod, 24 {degrees}C) and four FR timing treatments: no FR (Control), FR throughout production (FullFR), FR during early development only (EarlyFR), and FR during late development only (LateFR). LateFR increased marketable fresh weight relative to Control (244 vs 224 g) and reduced flowering incidence, whereas far-red supplied during early development reduced fresh weight (158 g) and increased flowering. The magnitude of the timing response differed among cultivars: switching from EarlyFR to LateFR recovered 0 % fresh weight in Amador but 107 % in Renegade and Rubino, with the largest penalties occurring in otherwise bolt-resistant cultivars. EarlyFR also increased total chlorophyll and reduced the chlorophyll a:b ratio. These results show that FR response in spinach is strongly conditioned by developmental stage and cultivar. Although LateFR received more total far-red than EarlyFR, it behaved like the Control, indicating that the penalty was set by far-red timing rather than dose. Treatment differences in bolting and yield tracked an estimated phytochrome photostationary-state deficit during early development: a phytochrome-deficit model markedly outperformed a cumulative-dose model ({Delta}AIC = 441), and the deficit x cultivar interaction was strong (p < 0.001), with bolt-resistant cultivars losing most yield when far-red coincided with the early developmental window. We therefore propose that FR should be treated as a genotype-dependent management variable rather than as a fixed spectral input, with late application and bolt-resistant cultivars offering the most favourable combination for vertical-farm spinach production. Framed within the breeders equation, the close match between the trial and production environment and the scope for shorter breeding cycles indoors suggest that genotype and far-red timing can be optimised jointly to accelerate genetic gain.

14
Bidirectional hybridization between Ulva prolifera and U. linza (Ulvophyceae, Chlorophyta): Evidence for compatibility and paternal chloroplast inheritance

Xu, Z. Z.; Zhao, J.; Jiang, P.

2026-07-22 genetics 10.64898/2026.07.18.739365 medRxiv
Top 0.2%
2.4%
Show abstract

Ulva prolifera and U. linza are closely related species, with abundant adult thalli and reproductive cells co-occurring extensively in time and space during the Yellow Sea green tides. Elucidating their hybridization compatibility is crucial for species delimitation, assessing interspecific gene flow, and evaluating the ecological impacts of green tides. Previous studies suggested asymmetric gamete compatibility (only U. prolifera mt+ x U. linza mt-), but lacked sex-linked markers to reliably identify hybrid diploids and their reproductive modes, and did not examine chloroplast inheritance. Here, we performed bidirectional crosses using sexual strains of different geographic origins from both parents, with sex-linked markers to quantify progeny genotypes, examine fertility and reproduction pathway of F1 hybrids, and trace chloroplast inheritance using the species-specific petB marker. Our results showed that: (1) F1 hybrids were obtained in both cross directions, with significantly higher frequency in the direct cross (U. linza mt+ x U. prolifera mt-), indicating no complete reproductive isolation in either direction; the biased compatibility likely reflected genetic background differences among the limited strain combinations in a single study. (2) A considerable number of germinated progeny arose from parthenogenesis of parental gametes. (3) F1 hybrids from both crosses could undergo meiosis to form gametes and develop into gametophytes; additionally, F1 from the reciprocal cross produced diploid spores for asexual reproduction, suggesting meiotic disturbance. (4) Chloroplasts were maternally inherited in selfing of U. prolifera parent, but in all F1 hybrids they were paternally inherited, indicating a potential reversal of the inheritance pattern due to interspecific hybridization. These findings provided new insights into the potential for genetic exchange between U. prolifera and U. linza.

15
VC1 and the production of vicine and convicine in the genus Vicia

Vottonen, L. L.; Chang, W.; Pöysä, M.; Lampi, A.-M.; Tanskanen, J.; Schulman, A. H.; Stoddard, F. L.

2026-07-14 plant biology 10.64898/2026.07.09.737524 medRxiv
Top 0.2%
2.4%
Show abstract

Many Vicia species contain vicine and convicine (VC), which limit the use of faba bean and some vetches in food and feed. The first step in VC biosynthesis in V. faba is shared with the riboflavin pathway and attributed to VC1, a member of the ribAB family. Since riboflavin is ubiquitous to life, we examined the distribution of VC production in genus Vicia. Three accessions of each of 33 Vicia species were grown in glasshouse conditions to provide fresh seeds for VC analysis and leaves for DNA analysis. PCR was used to amplify fragments of the VC1/ribAB gene for sequencing, and these sequences were used to create a phylogenetic tree. COX1 and ITS2 sequences were used for examining the nucleotide diversity in the subgenera. VC and DNA sequences consistent with VC1 were found only in members of subgenus Vicia. In V. lathyroides, VC1 was present but no VC was detected. There was less sequence diversity in VC1/ribAB sequences of subgenus Cracca than in those of subgenus Vicia, suggesting that ribAB remained under stricter purifying selection than VC1. VC1 is confirmed as a prerequisite for the presence of VC, and the gene and its products are restricted to subgenus Vicia. HighlightThe favism-causing factors of vetches and faba bean, vicine and convicine, depend on the presence of the VC1 variant of the ribAB gene, which is found in only one subgenus.

16
Sunrise and sunset times are the main factors that determine the flowering time of photoperiod-sensitive sorghum

Clerget, B.; Sidibe, M.; vom Brocke, K.; Raharinivo, V.; Ortiz, D.; Trouche, G.

2026-07-08 plant biology 10.64898/2026.06.12.731875 medRxiv
Top 0.2%
2.4%
Show abstract

Crop photoperiodism models assume that flowering time is primarily controlled by daylength, yet many field observations contradict this view. We previously proposed an alternative framework integrating daily changes in sunrise and sunset times (dSR and dSS). Variety trials in Madagascar and in Argentina supported this concept: mid-late sorghum varieties from the northern hemisphere flowered late or very late when sown in November and December, consistent with the higher dSR/dSS values of the southern hemisphere summer. One Malian variety, sown monthly over six years in West Africa, exhibited high interannual variability in flowering time when sown between November and February. This revealed that up to four photoperiodic responses -- two quantitative and two qualitative, occurring at different times of the year -- may coexist within a single late photoperiod sensitive variety. All responses use only dSR and dSS cues. The qualitative responses are triggered by an internal phasic coincidence, which is set by a linear relationship between dSR and dSS at the onset of plant photoperiod sensitivity, and between dSR+dSS at panicle initiation. The research model fitted data from 28 varieties grown in Mali well. It also accurately fitted the duration to PI observed in three varieties sown at tropical and temperate latitudes. HighlightThe seasonal photoperiodic adaptation of flowering time in sorghum plants may rely on several signal transduction pathways regulated by sunrise and sunset times rather than day length.

17
Evaluating crop models for future climate scenarios: wheat yield predictions using APSIM and STICS under combined CO2, warming, and water deficit conditions

Severini, A. D.; Gawinowski, M.; Bancal, M.-O.; Launay, M.; Deswarte, J.-C.; Chenu, K.

2026-06-10 plant biology 10.64898/2026.06.07.730737 medRxiv
Top 0.2%
2.4%
Show abstract

Crop models are essential for predicting climate change impacts on agriculture, yet their validation under multi-stress conditions remains limited. This study evaluated two widely-used wheat models, APSIM and STICS, using data from three Free-Air CO2 Enrichment (FACE) experiments (USA, Germany, Australia) combining elevated CO2 (eCO2), water deficit, and warming. Environmental characterisation using simulation-based stress indices revealed that intended "controls" frequently experienced hidden heat and water stress, meaning models were calibrated on crops already undergoing physiological adjustments. Evaluation of simulated yield and components revealed a clear hierarchy in prediction errors (RRMSE): unlimited conditions (3-9%) < single stress (4-27%, with a need to improve response to heat stress) < combined stress (17-123%). Elevated CO2 generally increased prediction uncertainty for crops experiencing water stress. Our results suggest that current stress functions from the models fail to capture the synergistic coupling between drought and heat stress. This highlights the urgent need for more mechanistic modelling to improve the reliability of climate change impact assessments.

18
Plant regulator of flower bud differentiation in in vitro plants of Cymbidium tortisepalum var. longibracteatum with TDZ as the key initiator

WU, C. q.; ZHAO, T.; CAO, F.; LI, H.; LI, J.

2026-07-30 plant biology 10.64898/2026.07.29.741633 medRxiv
Top 0.2%
2.4%
Show abstract

BackgroundCymbidium tortisepalum var. longibracteatum is a Class II nationally protected endangered plant in China with significant economic value. AimThis study aimed to establish an efficient in vitro flowering system and identify key hormonal factors regulating flower bud differentiation. MethodsOrthogonal experiments and hormone treatments were performed to evaluate the effects of TDZ, 6-BA, NAA, IBA, PP333, ABA, and GA3 on flower bud induction in tissue-cultured plantlets. ResultsTDZ was identified as the key initiator of flower bud differentiation, as no flower buds were induced in TDZ-free treatments. The optimal hormone combination was 1.0 mg{middle dot}L-1 6-BA + 0.4 mg{middle dot}L-1 TDZ + 0.8 mg{middle dot}L-1 NAA + 1.6 mg{middle dot}L-1 IBA, achieving a flower bud induction rate of 23.33% and a normal flower bud rate of 20.00%. PP333 inhibited flower bud differentiation but reduced malformation; at 0.2 mg{middle dot}L-1 PP333, the normal flower bud rate reached 20.00%. GA3 pretreatment resulted in a flower bud induction rate of 12.04%, whereas ABA pretreatment showed no significant promoting effect on flower bud induction. ConclusionsThis study provides an efficient in vitro flowering system and key hormonal parameters for shortening the breeding cycle and elucidating orchid flowering mechanisms.

19
Yield losses associated with peanut smut incidence in Argentina: a quantitative synthesis across field studies

Cazon, L. I.; Gonzalez, N. R.; Del Ponte, E. M.; Costa de Carvalho, A. C.; Asinari, F.; Camiletti, B. X.; Paredes, J. A.

2026-08-11 plant biology 10.64898/2026.08.11.744131 medRxiv
Top 0.2%
2.1%
Show abstract

Peanut smut, caused by Thecaphora frezzii, is an important constraint to peanut production in Argentina, but quantitative estimates of yield losses across environments remain limited. We quantified the relationship between disease incidence and kernel yield using 922 observations from 26 field studies conducted in Cordoba, Argentina, between 2021 and 2025. Study-specific incidence-yield relationships were analyzed using linear regression, random-effects meta-analysis, and linear mixed-effects models. Peanut smut incidence was consistently associated with yield reduction across studies. The estimated damage coefficient ranged from 24.2 to 28.7 kg ha-{superscript 1} per 1% increase in disease incidence, corresponding to a relative yield reduction of 0.74-0.87% of attainable yield. In contrast, attainable yield varied markedly among studies, ranging from 1,370 to 5,409 kg ha-{superscript 1}. Although an exploratory segmented analysis suggested a breakpoint near 12% incidence, subsequent moderator analyses, study- specific regressions, and normalized response curves provided no evidence of a biologically meaningful change in the damage coefficient across incidence or yield classes. These results indicate that differences among environments were primarily associated with attainable yield rather than with changes in the magnitude of disease-associated yield loss. The resulting damage function provides a quantitative basis for yield-loss assessment and disease management in peanut.

20
Small variations in temperature and photoperiod under field conditions in South America modified spring barley phenology and foliar development

Mastandrea, N. F.; Quero, G. E.; Castro, A. J.

2026-08-28 plant biology 10.64898/2026.08.27.747666 medRxiv
Top 0.2%
2.1%
Show abstract

Context. Barley production requires advanced knowledge of its response to changing environmental conditions in order to keep it competitive and sustainable. Aim. Advance in the understanding of barley phenology and foliar development under South American field conditions. Methods. 8 spring barley genotypes with differential phenology were studied in four field experiments under different temperature (through years and sowing dates) and photoperiod (through sowing dates) conditions. Time to anthesis, emergence to onset of stem elongation, stem elongation to anthesis, photoperiod response (PR) in these three traits, number of final leaves at anthesis (FLN) and phyllochron were measured. Key Results. Time to anthesis and its subphases were shorter in late plantings but under similar photoperiod, temperature increased them. Cultivars have differential responses but with magnitude interactions and not crossover ones. Cultivar effects defined PR with no interaction with year (temperature). Temperature and photoperiod affected FLN, phyllochron and their relationship with time to anthesis. Under the shorter photoperiod, FLN and phyllochron were negatively correlated, FLN was higher in the warmer year and positively correlated with time to anthesis while phyllochron was not affected by temperature and had no correlation with time to anthesis. Under longer photoperiod, phyllochron was higher in the warmer year and time to anthesis was positively correlated with both FLN and phyllochron. Conclusions. Cultivar basal thermal requirements and PR were consistent under the different studied conditions. Changes in temperature and photoperiod affected the relationship between time to anthesis, FLN and phyllochron suggesting that, although the three traits are arithmetically related, environmental conditions affect their balance.